Carbon Nanotube Thermal Transport: Ballistic to Diffusive 
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We propose to use lo/(lo + L) for the energy transmission covering both ballistic and diffusive 
regimes, where lo is mean free path and L is system length. This formula is applied to heat 
conduction in carbon nanotubes (CNTs). Calculations of thermal conduction show: (1) Thermal 
conductance at room temperature is proportional to the diameter of CNTs for single-walled CNTs 
(SWCNTs) and to the square of diameter for multi- walled CNTs (MWCNTs). (2) Interfaces play 
an important role in thermal conduction in CNTs due to the symmetry of CNTs vibrational modes. 
(3) When the phonon mean free path is comparable with the length L of CNTs in ballistic-diffusive 
regime, thermal conductivity k goes as L a . The effective exponent a is numerically found to decrease 
with increasing temperature and is insensitive to the diameter of SWCNTs for Umklapp scattering 
process. For short SWCNTs (< 0.1 /im) we find a ~ 0.8 at room temperature. These results are 
consistent with recent experimental findings. 
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Carbon nanotubes (CNTs) have been reported to have 
remarkable electrical, mechanical, and thermal proper- 
ties, making them the ideal material for various engineer- 
ing applications. 1,2 Of particular interest in this letter is 
thermal conduction of CNTs. For example, the anoma- 
lous enhancement of thermal conductivity has been ex- 
perimentally observed in CNTs composites^. However, 
many conflicting results of CNTs thermal conduction 
have been reported both theoretically and experimen- 
tally. The values of CNTs thermal conductivity rang- 
ing from 30W/mK to 6600 W/mK have been reported 
at room temperatureA^SiLSiS These drastically differ- 
ent values make it necessary to clarify which is reliable 
and what are the conditions for obtaining such values. 
Thermal conduction in CNTs may differ from the pre- 
dictions of Fourier's law based on bulk materials be- 
cause the phonon mean free path (MFP) can be com- 
parable to the length of CNTs. Ballistic phonon ther- 
mal transport behavior in CNTs has been observed in 
molecular dynamic (MD) simulations^ and recently in 
experiments—*^ However, so far as we know, no satis- 
factory theory has been able to cover both ballistic and 
diffusive region. Landauer formula— takes care of bal- 
listic thermal transport, while Boltzmann equationsi^*^ 
cannot go over to purely ballistic transport. 

In this letter, we propose a formula for thermal conduc- 
tion covering the range from ballistic to diffusive trans- 
port. With this formula, we discuss the CNTs thermal 
conduction's dependence on the tube diameter and tube 
length, as well as interface effect. 

For thermal transport in quasi-one-dimensional sys- 
tem, the thermal conductivity can be written asi£ 

n,v n >0 

where L and S are the length and cross-section area of 
the system, T n (q,u> n ) the energy transmission for the n- 
th branch wave at longitudinal momentum q, angular 
frequency uj n , and f(uj n ,T) is the Bose- Einstein distri- 



bution at temperature T. v n is group velocity along the 
direction of thermal transport given by v n — duj n /dq. 
The integration over momentum q is within the first 
Brillouin zone. The central problem is to calculate the 
energy transmission T n (q,uj n ). We propose that energy 
transmission can be calculated for three different ther- 
mal transport regimes according to the relation between 
the MFP l (q) and the length of the system L. (1) Bal- 
listic regime lo 3> L. The MFP is much larger than 
the length of the system. Thermal transport in this 
regime can be described by Landauer quantum trans- 
port formula. 1 — There are two equivalent methods for 
calculating energy transmission resulting from boundary, 
conjunction, or defect scattering from atomistic point of 
view: the scattering boundary mode match method^ 
and the non-equilibrium Green function method— (2) 
Ballistic-diffusive regime where lo ~ L. The MFP is 
comparable with the length of the system. In this regime, 
the energy transmission is approximated by T n (q,u> n ) — 
Iq/(L + Iq), as has been shown in Ref. H?] for electronic 
transport. (3) Diffusive regime Iq <C L. The MFP is 
much smaller than the length of the system. In this 
regime, T n (q,uj n ) ~ Iq/L. So thermal conductivity in 
this limit is given by the well-known Boltzmann-Pcicrls 
formula K = £ £/ ^hw n (q)v n ^(q). 

n 

Diameter dependence. Here we assume that the ther- 
mal transport in CNTs is in the ballistic regime. To 
simplify the discussion and catch the essential charac- 
ter of ballistic thermal transport in CNTs, we assume 
the energy transmission T n (q,u> n ) ~ 1. This assump- 
tion is justified for short CNTs with few imperfections 
at moderate temperature a 9 i 10 . Several different defini- 
tions for the cross section for CNTs have been used in 
Ref. |4 5 6 7 8 9, because it is not well defined. To cir- 
cumvent this ambiguity, we would like to use more exper- 
imentally oriented physical quantities: thermal conduc- 
tance Gth = kS/L to measure the property of thermal 
conduction of CNTs. In our computation, phonon dis- 
persion for single-walled CNTs (SWCNTs) is calculated 
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FIG. 1: Thermal Conductance for SWCNTs. (1A) Differ- 
ent chirality SWCNTs as a function of temperature. (IB) 
Different diameter SWCNTs at T = 50, 100, 200, 300 K re- 
spectively. 

using the method and force constants in Ref. Q. The 
group velocity is calculated from the phonon dispersion 
by the method in Ref. 0. The results of thermal con- 
duction for different chirality and diameter SWCNTs at 
different temperatures are illustrated in Fig. ^ It can 
seen from Fig. ^ that at low temperature, for example 
T = 50 K, thermal conductance almost does not change 
with diameter or chirality. At the extremely low tem- 
perature limit, all CNTs give the quantized universal 
thermal conductance^ 4-7r 2 fc 2 3 T/3/i, which corresponds 
to the four channels of quantum transport for acous- 
tic phonons. When temperature increases, more optical 
phonon channels will contribute to thermal transport. 
However, it is contrary to the intuitive conclusion^ that 
thermal conductance at high temperature T will be pro- 
portional to the total numbers -/V p h of phonon channels 
in CNTs. Thermal conductance at high temperature is 
proportional to the diameter d of CNTs, irrespective of 
their chirality. For example, the ratio of channel number 
SWCNTs (6,4) to (6,0) is / N^ 0) » 6.33, while 

their diameter ratio is d^ 6 ' 4 ' /(V- 6 ' ' ~ 1.45. It can be 
seen from Fig.^that the ratio of thermal conductance at 
300 K is about 1.45, not 6.33. This proportionality to di- 
ameter is implicitly demonstrated in Ref. |1J by counting 
the number of phonon branches at frequency to to calcu- 
lating the transmission. This simple result can be easily 
understood if we approximate the integrand in Eq. 
by a constant, obtaining the result that Gth is propor- 
tional to the total number of modes divided by the unit 
cell length, which is proportional to d. For MWCNTs, if 
thermal conductance above certain temperature can be 
added from the contribution of each shell independently, 
the total thermal conductance will be oc J rdr oc d 2 . The 
dissipation power results in Ref. [fci have a tendency of 
d 2 . It should be noted that the above results only charac- 
terize the thermal conductance in pure ballistic regime. 
In actual experiments, the situation is more complicated 
for thermal conductance in sub-micrometer lengths. For 
example, the contact interface and phonon scattering are 
inevitable in experiment. We will discuss these effects on 
thermal conduction in CNTs. 

Interface effect. Interface effect on CNTs thermal 
transport has been studied through effective medium 



theory^ and empirical random walk simulation^ 19 How- 
ever they are only rough models, not from first principles. 
Here we consider the interface effect from a lattice point 
of view. We consider an extreme case, where the semi- 
conductor nanotube junction structure (11,0) and (8,0) 
is first constructed by a geometrical method as in Ref. 1. 
The structure is optimized by a second-generation Bren- 
ner potential 2 - to let the atoms get their equilibrium 
positions. The force constants are derived numerically 
from the same potential. The phonon dispersion is calcu- 
lated from these linearized force constants. Four acoustic 
branches are considered for energy transport: the longi- 
tudinal mode (LA), doubly degenerate transverse mode 
(TA), and the unique twist mode (TW) in nanotubes. 
Following the scattering boundary method^ we calcu- 
lated the energy transmission across the CNTs conjunc- 
tion. The details of calculation are presented in Ref. fl5L 
The reflected and the transmitted waves across the junc- 
tion for the incident LA mode waves are both only LA 
modes. We think that this is due to the high symmetrical 
properties of atomic motion for nanotubes. The LA and 
TA modes are common symmetrical motions for both the 
left and right lead. So they propagate through the con- 
junction. In contrast, the transmissions of the TW mode 
and many other optical modes are nearly zero or very 
small. From these results, we consider the interface be- 
tween CNTs and substrate will influence the experimen- 
tal measured CNTs thermal conductivity greatly which 
has been indicated in experimental result a 8 i 10 . We argue 
that chemical functional reorganization of the CNTs ends 
or proper choice of the more symmetry-shared substrate 
matrix will improve thermal transport of CNTs. 

Length effect. Purely ballistic thermal transport will 
make thermal conductivity k diverge linearly with the 
length L of CNTs. Anomalous thermal transport for ID 
nonlinear lattice predicts^ that k will diverge with L as 
k oc L a . For ID nonlinear lattice with transverse motion, 
which is similar to CNTs, it is found a = 1/3 through 
MD simulation and mode-coupling theory^ 2 . A few re- 
sults of tube length effect on thermal transport in CNTs 
have been reported recently through MD simulation 7 - 9 . 
A length dependence of the thermal transport in CNTs 
is also observed in experiments^. It can be seen from 
Eq.Q that thermal conductivity will depend on the 
length of CNTs in ballistic-diffusive regime. 

The problem to calculate thermal conductivity from 
Eq. m in ballistic-diffusive regime is to compute MFP, 
which depends on frequency and wave length of the vi- 
bration modes, as well as temperature. Assuming that 
the scattering mechanisms are not coupled, MFP can 
be written 2 ^ as 1/Iq = ^jl/Zi, where the index i de- 
notes various interaction phonon processes. In CNTs, 
the phonon scattering process may be different from that 
of bulk and is not well known. Here, for simplicity, we 
first consider the Umklapp scattering. Umklapp phonon 
scattering Ifj 1 in CNTs has been shown linear T depen- 
dence at high temperatures through experiment 8 and is 
usually given a lo 2 dependency on the frequency^* So 
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we take the form Ijj = A/(uj 2 T) for Umklapp scattering 
in CNTs, where A is taken to be the same as derived 
for graphene in Ref. |3 The phonon dispersion is cal- 
culated as in Ref. Q. The cross-section area S is denned 
as S = 7rcP/4. Results are illustrated in Fig. |3] through 
Eq. JJJ for ballistic-diffusive regime in CNTs. We should 
note that the exponent a in Fig. GIC,D) is an effective 
exponent in a limited length region, not the asymptotic 
one> 2i i 22 j^. can kg geen f rom pig 2A that thermal con- 
ductivity shows different length dependence scaling L a 
for short CNTs (< 0.1/xm) and long CNTs (> 10/zm), 
denoted as as and ai respectively, when Umklapp scat- 
tering in CNTs is considered in ballistic-diffusive regime. 
This is due to the relation between MFP and the length 
of CNTs. Average MFP l (T) = f l Q {uj,T)D{uj)duj is 
calculated from the phonon dispersion relation for CNT 
(6,6), where D(lo) is normalized density of states, satis- 
fying J D(u>)du) = 1. Figure [33 indicates that average 
MFP decreases rapidly from 0.25 /mi to 0.06 /tm, when 
temperature goes up from 50 K to 200 K. We also calcu- 
lated average MFP for other chirality CNTs and found 
little difference from that of CNT (6,6). Fig. 2C shows 
the change of as and cti with temperature. For exam- 
ple, when T = 50 K, the average MFP is as large as 
0.25 /jm and so the as is almost equal to f. When tem- 
perature increases, the MFP becomes shorter, so as and 
ai both decrease with temperature. The dependence of 
as and ql on CNTs diameter is illustrated in Fig. 2D. 
as almost does not change with the chirality of CNTs 
and gives as ~ 0.83 at room temperature, while a^ de- 
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creases slowly with increase of the diameter of CNTs, and 
gives its value about ~ 0.35 for CNTs with diame- 
ter larger than Iran at room temperature. This kind of 
temperature dependence and chirality dependence agrees 
with MD simulation resultsi^ It is worth noting that, 
for CNTs with length less than 0.1 fim, a = 0.32 and 
a = 0.40 for CNT (5,5) at room temperature are reported 
by non-equilibrium MD simulation^, while a ps I is ob- 
served using equilibrium MD simulation^. For MWCNT 
with length less than 0.5 /xm, roughly a ~ 1 is reported^ 
We think that the difference for non-equilibrium MD and 
equilibrium MD comes from the boundary condition be- 
cause in short CNTs the result should be sensitive to the 
boundaries. 

Next we discuss other phonon scattering mechanisms. 
If the frequency dependence of MFP is assumed as 1/Zo oc 
uj t , then at high temperature the length dependence of 
thermal conductivity given by Eq. Q can be estimated 

as k ~ J 1 i lo ^_ 1 i L duj oc L~ ~, that is a = (r — l)/r. 
For Umklapp scattering process, r — 2 and a = 0.5, 
which corresponds to the effective exponent at length 
about 1/xm. For defect scattering mechanism, we have 
1/Zo oc ut 4 , a — 0.75. Thus, it can be argued that 
strong frequency dependence scattering (large r) will not 
have much help in eliminating thermal conductivity di- 
vergence, while weak frequency dependence scattering 
will contribute to eliminating this divergence. Finally, in 
the very long tube limit, n approaches a constant within 
our theory. 

Conclusions. Crossover in thermal conduction from 
ballistic to diffusive is illustrated. Thermal conductance 
at room temperature is found oc d for SWCNTs and 
oc d 2 for MWCNTs. Interfaces plays an important role 
in thermal conduction in CNTs due to the high sym- 
metrical property of CNTs vibration modes. Possible 
ways to improve thermal conduction in CNTs are sug- 
gested. In ballistic-diffusive regime, thermal conductiv- 
ity k behaves as L a . The effective exponent a is nu- 
merically found decreasing with increasing temperature 
and insensitive to the diameter of SWCNTs for Umk- 
lapp scattering process. The possible mechanism for the 
reduction of divergence for thermal conductivity is also 
discussed. Although the formula is still phenomenolog- 
ical, it does cover ballistic and diffusive regimes with a 
smooth crossover and gives a reasonably simple picture in 
the whole temperature range. This work is supported in 
part by a Faculty Research Grant of National University 
of Singapore. 



FIG. 2: Length divergence in ballistic-diffusive regimes: 
(2A) Thermal conductivity for different chiralities SWC- 
NTs at 300 K; (2B) Average mean free path length for nan- 
otube (6,6) at different temperatures; (2C) Change of power 
law exponent for short and long CNTs with temperature; 
(2D)Change of power law exponent for short and long CNTs 
with diameter; 
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